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ABSTRACT 

From the launch of the Fermi Gamma-ray Space Telescope to July 9, 2010, the Gamma-ray Burst 
Monitor (GBM) has detected 497 probable GRB events. Twenty-two of these satisfy the simultaneous 
requirements of an estimated burst direction within 52° of the Fermi Large Area Telescope (LAT) 
boresight and a low energy fluence exceeding 5 /ierg/cm 2 . Using matched filter techniques, the spatially 
correlated Fermi /LAT photon data above 100 MeV have been examined for evidence of bursts that 
have so far evaded detection at these energies. High energy emission is detected with great confidence 
for one event, GRB 090228A. Since the LAT has significantly better angular resolution than the 
GBM, real-time application of these methods could open the door to optical identification and richer 
characterization of a larger fraction of the relatively rare GRBs that include high energy emission. 
Subject headings: Gamma-ray burst: general 



1. INTRODUCTION 

One of the more surprising results of the Compton 
Gamma-Ray Observatory (CGRO) was the EGRET dis- 
covery of an 18 GeV photon associated with GRB 940217 
(Hurley 1994). About a half-dozen bursts were seen over 
the course of the CGRO mission with photons above 100 
MeV (Catelli 1998, Dingus 2003). Since the GRB spec- 
tral energy distribution at lower energies has been well 
characterized by a modified power law with peak fluxes at 
energies of the order of 200 KeV, the existence of photons 
at energies 10 times higher puts a significant constraint 
on any viable model of the GRB phenomenon. This has 
been a subject of great interest for missions that followed 
EGRET. Prior to launch of the Fermi Gamma-ray Space 
Telescope, it was possible to speculate that the LAT in- 
strument would detect more than 200 GRB events per 
year (Dingus, 2003). In the two year period since the 
launch of the Fermi Gamma-ray Space Telescope, the 
Gamma-ray Burst Monitor (GBM) has reported approx- 
imately 475 GRBs, ie. a rate of about 250 per year. 
Over essentially the same period, only 17 bursts have 
been identified by the Fermi/LAT. We now see that the 
range of GRB photon energies extends over a scale of 
10 6 but the physical dynamics of these phenomena are 
still not understood. This is coupled to the question 
of whether high energy photons are associated with all 
GRBs or only with a small sub-class. Since the Fermi 
mission is unlikely to be duplicated any time soon, there 
is some urgency to assuring that the maximum informa- 
tion is being extracted from this valuable facility. Thus, 
our group has set about developing techniques for en- 
larging the number of gamma-ray bursts identified with 
high energy photon emission, ie. above 100 MeV. 

The first result of this effort has established the 
correlation of two Swi ft /XKT -localized bursts, GRB 
080905 A and GRB 091208B, with high energy photons 
in the Fermi/LAT detector (Akerlof el al. (2010), here- 
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after A10). The statistical technique employed is the 
matched filter method, most familiar to those detect- 
ing signals in the time domain. The underlying assump- 
tion is that the characteristics of both the signal and 
background are a priori known functions of one or more 
variables. Since the matched filter maximizes the signal- 
to-noise ratio, moderate departures from optimality de- 
grade the filter performance relatively slowly, making this 
a valuable tool for investigating the possible existence of 
faint signals. The details of the filter algorithm are ex- 
plicitly described in A10. In this paper, we take the next 
harder step of dropping our reliance on precision burst 
coordinates provided by Swift or other similar high res- 
olution instruments. Instead, we use the approximate 
localization of the Fermi /GBM to map a region of in- 
terest on the Fermi/LAT field of view. By identify- 
ing high energy photon clusters, provisional burst coor- 
dinates can be determined with significantly smaller er- 
rors than available from the GBM. From there, the burst 
identification follows along lines set out in A10. 

2. SAMPLE SELECTION 

As a first step in this program, a list of all GBM trig- 
gers was obtained from the fermigbrst catalog main- 
tained by the Fermi Science Support Center 3 . The cat- 
alog contains 497 GRB triggers from launch to July 9, 
2010. This list was cross-matched with Table 1 in Guetta 
& Pian (2009) and Table 2 in Guetta et al. (2010) to 
identify the burst GCN designations and the low energy 
fluences. For triggers occuring after February 18, 2010, 
fluences were obtained from individual GCN circulars. 
GBM triggers were also checked against XRT locations 
from Swift 4 to remove events already considered in A10. 

Using data from the Fermi spacecraft attitude file, we 
further selected those triggers with a boresight angle less 
than 52° and an estimated GBM error circle less than 
10°. Events without GBM fluence information or previ- 
ously claimed LAT detections 5 were also discarded. Ap- 

3 http:/ /heasarc. gsfc.nasa.gov/W3Browse/fermi/fermigbrst. html 

4 http://heasarc.gsfc.nasa.gov/docs/swift/archive/grb_table/ 
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TABLE 1 
List of 22 GBM trigger GRBs 



GRB 


Trigger 


RA 


Dec 


Sqbm 






(°) 


(°) 


/ierg/cm 2 



9.2 

5.0 
10.9 
19.4 

9.6 
12.0 

5.2 
22.3 

9.0 

6.1 

7.5 
11.4 

8.1 
30.0 
102.0 

6.4 
11.4 
30.2 
10.0 

7.7 
12.3 

7.1 



plying a final cut on GBM fluence (8 - 1000 KeV) requir- 
ing greater than 5.0 /xerg/cm 2 reduced the number to 22 
events (see Table 1). These are termed the "GBM" data. 
464 additional fields were taken at random on the sky 
with similar criteria to study the background behavior 
and are identified as the "random" data in the following 
text. 

3. SIGNAL DETECTION TECHNIQUE 

The core task of this search procedure is the identifi- 
cation of triplet clusters of photons in the Fermi /LAT 
instrument whose spatial accuracy is considerably better 
than the GBM. The set of photon data for each candi- 
date burst is confined to lie within a 16° cone angle of the 
GBM direction and a time window extending from zero 
to 47.5 s after the GBM burst trigger. The procedure 
first computes a signal weight for each photon pair based 
on photon energy, detection time, photon event class and 
angular separation relative to the expected LAT PSF er- 
rors. The photon pair weights are subject to a weak 
threshold cut designed to avoid combinatorial overload 
should large photon numbers be encountered. In prac- 
tice, this was not a severe problem and can be ignored. 
The formula for the pair weights is given by: 

Qij = Wi ■ Wj ■ A^, (1) 

where 

Wi = w E (i) ■ w t (i) ■ w c (i) ■ 4:TT<jp SF (Ei) 7 (2) 

Sii ~^^PSFm+^PSFi.EiV (4) 

and Qij is the angle between the z'th and j'th photon. 
The definitions of we, w t and w c can be found in equa- 
tion 1,3 and 4 of A10. 

The next step is to link photon pairs so that the three 
pairs, {i, j}, {j, k} and {i, k}, become identified as the 



triplet, {i, j, k}. The triplet weight value is computed by 
the formula: 

Rijk = (uii ■ wj ■ w k ■ ■ A jk ■ A 4fc )3 (5) 

The triplet weights are ranked by value and the set is 
pruned by the condition that a triplet element, Rum, is 
removed if Ru m < Rijk ■ This leaves a set of triplet clus- 
ters, each with a discrete complement of three photons. 
For each of these clusters, a PSF-weighted estimate of 
the burst direction is performed and the matched filter 
angle weight, wq is computed with respect to this vec- 
tor as described in equation 2 of A10. At this point, 
an event weight for each cluster is computed by the for- 
mulas given in A10 with one small modification. In the 
scheme described here, the GRB direction is not initially 
defined with any precision. Thus, it is inappropriate to 
include a l/cr PSF factor for all values of we. For the 
highest energy photon in each triplet, the A-K(jp SF fac- 
tor is removed to reflect that this leading photon plays 
the principal role in fixing the apparent GRB direction. 
Although the calculations carry each cluster through the 
same computational path, the expectation is that the 
cluster with the highest matched filter weight is the most 
probable identification. 

The 22 GBM fields described earlier were the target 
of this investigation. We recognized that the most con- 
vincing argument for a true LAT identification should 
rely on the statistical distributions for the matched filter 
weights in LAT fields with similar characteristics. To in- 
crease that number as much as possible, the LAT fields of 
view were segmented into 12 circular tiles embedded on a 
spherical surface. Each tile subtends a cone with a half- 
angle of 16.0°. This tiling scheme was applied to both the 
GBM and random field data sets to realize 182 and 3440 
independent directions in space satisfying all the criteria 
described previously. Taking advantage of the fact that 
each field observation was blocked into a 250-s segment, 
the number of independent observations was multiplied 
by five by regarding each 50-s time slice as a separate 
sample. Thus, there are 910 background measurements 
taken from LAT observations obtained simultaneously 
with the candidate GBM fields and an additional 17200 
samples taken under similar but not identical conditions. 
One particular concern for an analysis of this type is that 
false positives will selectively occur as the sample photon 
rate rises substantially above the mean. Evidence that 
this is not the case here is shown in Figure 1 which plots 
the cumulative distribution of the total number of pho- 
tons within the LAT field of view over a 250-s interval. 
These rates explicitly exclude contamination from pho- 
tons beyond the 105° zenith angle cut. As shown in the 
plot, the prominent GBM trigger event reported here is 
not associated specifically with fields with high ambient 
background rates. The similarity of the distributions for 
GBM and random fields also shows that the GBM data 
are not pathological as far as rates are concerned. 

4. RESULTS 

Our statistical localization and weighting scheme iden- 
tified one outstanding candidate for high energy photon 
emission, GRB 090228A. The best estimate for the prob- 
ability of such an occurence by chance alone was ob- 
tained by performing identical searches on random LAT 
fields with the same criteria. Thus, 11 out of 17200 ran- 
dom fields generated matched filter weights exceeding the 
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Fig. 1. — Cumulative distributions of LAT photon rates over 
250-s intervals for the GBM (blue) and random background (red) 
fields. These rates reflect the entire LAT FoV except for photons 
that lie outside the 105° zenith angle cut. The rate corresponding 
to the most prominent GBM event is indicated by the arrow. 

value for our candidate event. Multiplying by a trials fac- 
tor of 22 for the number of GBM localized fields consid- 
ered yields a false positive probability of 1.4%. To check 
that these correlations were simply not due to preferen- 
tially higher background rates for the GBM exposures, 
we also performed similar calculations for the LAT data 
confined to an average of 8 uncorrelated directions per 
exposure and five independent time intervals from the 
same GBM data sets. In this case, 2 fields out of 910 
exceed the candidate signal for a false positive rate of 
4.8%. The cumulative distributions are plotted in Fig- 
ure 2. The statistical similarity of the GBM off-axis and 
random field data demonstrates that the GBM data set 
is not correlated with anomalous environmental condi- 
tions such as higher cosmic ray background rates. A list 
of photons associated with this burst is provided in Table 
2. 

The GBM data for GRB 090228A is described in GCN 
8918 (von Kienlin, et al., 2009). According to this note, 
the burst was localized to a 1 — a accuracy of better 
than 1° with an additional systematic uncertainty of the 
order of 2.5°. The coordinate values obtained by the 
GBM group and this analysis are listed in Table 2. The 
GBM value generated some concern since our cluster 
finder position disagreed by 8.7°. Fortuitously as this 
manuscript was being drafted, a paper was posted to 
astro-ph (Guiriec et al., 2010) providing a burst direc- 
tion with an estimated accuracy of 0.2° and lying 0.5° 
from our own estimate. We believe that this establishes 
the validity of our identification to near certainty. In 
addition, the positive GRB correlation of Event Class 2 
and 3 photon rates discussed in A10 was observed for the 
ensemble of GBM fields as well. The photon clustering 
is easily observed in the sky map shown in Figure 3. 

5. DISCUSSION 

The one event identified in this paper establishes the 
validity of our statistical techniques to a level of near 
certainty. By using GBM triggers to guide the discovery 
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Fig. 2. — Complements of the cumulative distributions for C"l2 w i 
for 22 GBM fields (blue), 910 random fields obtained nearly si- 
multaneously with the GBM data (red) and 17200 random fields 
obtained at random times (green). 



TABLE 2 

GRB 090228A high energy photon list and celestial 

COORDINATE ESTIMATES 
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2.007 


1.692 


125.241 
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52.814 
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3.752 
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25.141 
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0.002 
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GBM a 


106.80 


-24.30 


^ 3.0 






LAT 6 


98.56 


-28.86 


0.18 


8.66 




IPN C 


98.30 


-28.40 


0.02 


0.51 



indicates diminished we for highest energy photon 
a von Kienlin et al. (2009) 
b this paper 
c Guiriec et al. (2010) 

d di-\ i is the angle between the spatial directions for the 
GBM and LAT or the LAT and IPN 

of photon clusters in the LAT, the phase space for find- 
ing counterparts can be reduced from hundreds of square 
degrees to one square degree or less. This makes a very 
significant difference for those seeking to identify GRB 
optical counterparts. If the algorithms used here could 
be adapted to the real time environment, the number of 
bursts with high energy associations could be increased 
appreciably. The additional computational load is negli- 
gible - about 30 ms per day. For such real time applica- 
tions, the high selectivity employed here is overkill - any 
identification that can be corroborated optically will suf- 
fice. Thus, effective signal-to-noise rates of the order of 
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unity are extremely valuable. As shown here, these tech- 
niques greatly enhance the dynamic range over which 
high energy radiation can be explored. It is not too out- 
rageous to claim that this is the equivalent of making the 
LAT three to ten times larger in size. 

As was noted in A10, the most surprising aspect of 
our recent work is the very small number of GRBs that 
can be positively identified with high energy emission 
despite the substantially lower fluence thresholds. This 
is a mystery which deserves serious consideration. We 
hope that raising new questions is sometimes more useful 
than answering old problems. 



LAT photon sky map for GRB 090228A 

Fig. 3.— Sky map of > 100 MeV photons for GRB 090228A. The 
diameter of each dot is proportional to its statistical weight. Thus, 
the largest diameters represent Event Class 3, etc. The dotted 
circles around each point indicate the 1 — a errors. The figure is 
centered on the nominal coordinates furnished by the GBM; the 
blue dot on the lower left shows the GRB coordinates computed by 
the cluster algorithm described in the text. The large green circle 
depicts the boundaries of the 16.0° cone that defines the fiducial 
boundaries for the cluster search. North is up and East is to the 
right. 
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